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ABSTRACT 

A  theoretical  and  experimental  investigation  has  been  performed 
on  the  determination  of  the  characteristics  of  two  impinging  plane 
turbulent  jets  for  the  purpose  of  understanding  the  performance  of 
proportional  fluid  amplifiers  and  momentum -exchange  devices.   The 
jet  deflection  angle  was  determined  through  the  application  of  the 
free-streamline  theory  for  two  normally  impinging  jets  with  arbitrary 
throat  widths  and  wall  set-backs.   The  deflection  angle  and  the  veloc- 
ity and  turbulence  distributions  across  and  at  various  distances  along 
the  jet  were  determined  through  the  use  of  a  hot-wire  anemometer  system. 
The  results  are  presented  in  terms  of  normalized  parameters  suitable  for 
comparison  with  the  theoretical  predictions.   The  deflection  angles 
predicted  theoretically  were  found  to  be  in  good  agreement  with  those 
obtained  experimentally. 
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I.   INTRODUCTION 

In  recent  years  considerable  interest  has  been  shown  in  the  research, 
development,  and  design  of  fluid  amplifiers  with  no  moving  parts.   These 
devices  are  finding  many  applications  in  various  fields  of  engineering, 
especially  in  control  systems.   The  main  feature  of  most  of  these  devices 
is  the  use  of  a  controlled  fluid  jet.   The  power  jet  or  the  primary  jet 
may  or  may  not  attach  to  one  or  the  other  sidewall  depending  on  whether 
the  amplifier  is  a  wall-attachment  type  bistable  amplifier  or  a  propor- 
tional amplifier.   In  the  former  type,  the  deflection  of  the  jet  from 
one  stable  to  another  stable  position  is  accomplished  by  a  combination 
of  momentum  exchange  and  pressure  build-up  in  the  separation  bubble. 
In  the  proportional  or  momentum-exchange  devices,  however,  the  deflection 
of  the  jet  is  accomplished  by  and  proportional  to  the  exchange  of  momen- 
tum between  the  primary  and  control  jets. 

Despite  the  simplicity  of  their  design,  the  complex  combinations  of 
fluid  dynamic  phenomena  which  take  place  inside  these  devices  make  it  very 
difficult  to  produce  a  mathematical  model  able  to  predict  their  opera- 
tion.  The  need  for  gaining  a  better  understanding  of  the  basic  phenomena 
involved  in  the  fluid  amplifiers  has  brought  a  new  interest  in  the  area 
of  small  flows  in  closed  passages  and  chambers,  and  many  new  works  on 
these  subjects  have  appeared  in  the  last  few  years. 

Some  work  on  the  research  and  development  of  fluid  amplifiers  and 
rate  sensors  has  been  undertaken  by  the  Department  of  Mechanical  Engineer- 
ing of  the  Naval  Postgraduate  School  through  a  research  grant  from  the 
Harry  Diamond  Laboratories  of  the  U.  S.  Army  Materiel  Command.   The  pre- 
sent work  is  part  of  this  general  research  program  and  deals  with  only  the 
proportional  amplifiers. 
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A  schematic  diagram  of  a  typical  proportional  amplifier  is  shown 
in  Fig.  I.   The  interaction  region  of  the  amplifier  is  shaped  to  pre- 
vent the  attachment  of  the  supply  (power)  flow  to  an  output  via  the 
Coanda  effect.   With  no  input  (control)  flow  present,  the  supply  flow 
strikes  the  splitter  and  divides  equally  through  the  two  outputs.   If 
pressure  is  applied  to  one  of  the  inputs,  the  resultant  control  flow 
strikes  the  supply  flow  which  deflects  more  of  the  supply  flow  to  the 
opposite  output.   The  deflection  of  the  supply  flow  varies  proportionally 
with  the  input  flow.   The  greater  the  input  pressure,  the  greater  is  the 
deflection  of  the  supply  flow  to  the  opposite  output.   If  the  input  pres- 
sure differential  is  reduced,  the  deflection  of  the  supply  flow  is  re- 
duced proportionately.   Since  a  variation  in  the  control  flow  conditions 
causes  a  larger  variation  at  the  receivers,  the  device  can  be  considered 
as  an  amplifier. 

One  can  introduce  parameters  such  as  differential  flow  gain,  a  dif- 
ferential pressure  gain,  and  differential  power  gain  to  describe  the 
static  characteristics  of  the  amplifier.   Until  recently,  these  char- 
acteristics were  determined  purely  experimentally  without  a  serious  ef- 
fort to  theoretically  predict  the  jet  deflection,  momentum  interaction, 
velocity  and  turbulence  distributions,  and  the  division  of  flow  between 
the  two  outputs  and  vents. 

In  an  attempt  to  lay  a  theoretical  foundation  for  the  prediction  of 
the  characteristics  of  proportional  amplifiers  with  one  or  two  control 
flows,  the  present  study  was  undertaken.   Furthermore,  in  order  to  partly 
simplify  the  problem,  two  two-dimensional,  normally-impinging  jets  were 
considered  with  different  jet  widths  and  side-wall  set-backs.   Obviously, 
such  a  configuration  (Fig.  2)  corresponds  to  a  simple  momentum  exchange 
proportional  amplifier  as  encountered  in  several  logic  devices  such  as  a 
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FIGURE  2.   SCHEMATIC  DRAWING  OF  THE  JET  ASSEMBLY 
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"half  adder"  and  "transverse  impact  Modulators"  [1],   The  understanding 
of  the  fluid  mechanics  of  such  a  jet  configuration  could  then  be  ex- 
tended through  the  use  of  similar  but  perhaps  more  complicated  theoreti- 
cal and  experimental  means  to  the  study  of  the  complete  proportional 
amplifier  shown  in  Fig.  1. 

In  connection  with  this  and  other  current  and  anticipated  investi- 
gations in  the  laboratory,  it  was  decided  to  study  the  impingement  of 
two  semi-confined  jets  between  two  parallel  plates,  with  special  regard 
to  the  influences  of  the  velocity  ratios  of  the  two  jets  and  the  wall 
set-back.   A  large  scale  model  was  used  having  dimensions  that  are  large 
compared  with  the  power  of  resolution  of  the  flow  measurement  instru- 
ments.  Velocity  and  turbulence  measurements  were  carried  out  through 
the  use  of  a  hot-wire  anemometer  system  for  various  velocity  combina- 
tions and  at  various  sections  along  the  jet.   The  theoretical  analysis, 
based  on  the  free-streamline  theory,  was  used  to  calculate  the  jet 
deflection  angle  as  a  function  of  the  velocity  ratio  of  the  jets  assuming 
the  fluid  to  be  inviscid. 
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II.   THEORETICAL  ANALYSIS  OF  THE  JET  DEFLECTION 

Many  plane  irrotational  flow  patterns  which  include  the  formation 
or  deflection  of  free  or  partly  bounded  jets  can  be  analyzed  satis- 
factorily by  the  theory  of  free-streamlines  by  Helmholtz  [2]  and 
Kirchhoff  [3].   The  results  obtained  may  be  significant  not  only  for 
the  corresponding  two-dimensional  flows  of  real  compressible  and  in- 
compressible fluids  but,  in  many  instances,  for  their  three-dimensional 
counterparts  as  well.   The  analysis  consists  of  the  definition  of  suc- 
cessive conformal  transformations  involving  a  hodograph,  or  velocity 
plane,  and  the  application  of  the  Schwartz-Christoffel  transformation 
[4] .   The  fact  should  be  emphasized  that  from  a  mathematical  point  of 
view,  the  direct  calculation  of  free-streamline  flows  is  restricted  to 
two-dimensional  irrotational  flows  of  compressible  and  incompressible 
fluids  which  are  free  from  gravitational  effects.   Hence,  the  forma- 
tion of  a  jet  by  two  two-dimensional  streams  impinging  obliquely  is 
well  suited  to  free-streamline  analysis  because  of  the  dominance  of 
inertia  and  pressure  in  the  establishment  of  the  resulting  flow  pat- 
tern and  because  of  the  negligible  effect  of  gravity  on  gas  jets.   As 
will  be  discussed  later,  however,  difference  in  the  free-streamline 
velocities  of  the  two  jets  enhances  the  shear  in  the  mixing  region  of 
the  two  jets.   The  larger  is  the  velocity  difference,  the  larger  the 
energy  dissipation  becomes.   The  free -strema line  analysis  may  be  ap- 
plied, however,  only  when  the  velocities  along  both  sides  of  the  jet 
are  assumed  to  be  identical  and  equal  to  an  unknown  free-streamline 
velocity.   Thus  the  region  in  which  the  momentum-exchange  and  the 
turbulent  mixing  take  place  does  not,  strictly  speaking,  follow  the 
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assumptions  made  in  the  free  streamline  theory.   Hence  the  validity 
of  the  analysis  depends  on  the  experimental  demonstration  of  the 
uniformity  of  the  velocity  profiles  at  the  formation  region  of  the 
resulting  let. 

For  the  development  of  the  following  analysis,  the  free-streamline 
velocities  are  assumed  to  be  identical  and  the  validity  of  the  assumption 
is  discussed  later  in  connection  with  the  comparison  of  the  theoretical 
results  with  those  obtained  experimentally.   Mindful  of  the  restrictions 
and  limitations  of  the  applicability  of  this  method,  the  solution  of  the 
problem  under  consideration  proceeds  as  follows: 

If  7   and  w  are  used  for  the  complex  variable  and  the  complex 
potential,  then  the  variable   £  defined  as  dw/dz  is  given  by 

C  =  -  u  +  iv  -  qe  (1) 

in  which  u  and  v  are  the  velocity  components  in  the  x  and  y  directions, 
o  is  the  magnitude  of  the  velocity  vector,  and  0  is  the  angle# of  its 
inclination. 

If  the  variable  Q     is  introduced, 

ur        a 
0  =  Ln  (  -  -y-  )  =    -   Ln  ~  +  i©  (2) 


the  flow  region  in  the  physical  plane  (Fig.  3a)  can  be  mapped  onto  a 
corresponding  region  in  the  d   -plane,  as  indicated  in  Fig.  3b.   Either 
C   or  0  must  be  related  to  w  by  means  of  known  methods  through  still 
another  complex  variable  designated  as  t.   This  relationship  completes 
the  rather  complicated  functional  analysis  because,  as  has  already  been 
seen,  either  C  or  0  can  be  expressed  in  terms  of  dw/dz. 

Any  polygonal  boundary  in  one  complex  plane  can  be  transformed  into 
the  entire  real  axis  or  another  through  application  of  S  hwartz-Christof fel 
theorem,  the  interior  of  the  polygon  being  transformed  into  either  the 
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upper  or  the  lower  half  of  the  second  plane  [4].   Thus  the  outline  in 
the  CI   -plane  is  directly  transformable  into  the  real  axis  of  the  t-plane 
Therefore,  the    fi-plane  shown  in  Fig.  3b  may  be  transformed  to  the  t- 
plane  shown  in  Fig.  3c.   The  transforming  function  can  be  written  as 
follows: 

n  =    •         mt  (3) 

J  /(t+i)  (t-h) 

Upon    integration,    one   has 


CI  =    2MLn  C   (vt+1    +  vt-h)  (4) 

in  which  M  and  C  are  found  to  be  1/2  and  1/  I  +  h  by  substituting  into 
the  equation  the  coordinates  of  the  points  B  and  D.   The  final  form 
of  Q  thus  reduces  to 


Q  =    Ln  v't+I  +  /t^h  (5) 

v'h+1 

In  particular,    it  will   be  noticed   that  at    t  »  0,    fl    ■  ijJ         (at  point 

C)   and  one  has 

.  a         T        1    +    iyfr  .     .        -1  rr- 

ip  =   Ln =    i    tan     vh 

v'l+h 
or 


tan    /9  =  ■Jh 


(6) 


Obviously,    the   angle   of   the  deflection  of   the  main   jet,    Oi    ,    is   related 
to  $  by    the    simple    relationship 

a  +   j3  =  ff/2  (7) 
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FIGURE  3.   DEFINITION  OF  VARIOUS  FLOW  PLANES 
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The  velocities  are  related  to  various  parameters  in  the  t-plane 
via  the  Q     transformation.   At  A,  for  example, 


u  =  Ln  rr-   =  Ln  * 


Ua         ^/h~+l 


or 


Ua    ,       Jh+1  (8) 


uj        -vlc+1  +  ^/k^n 


Likewise,  at  E 


0  =  Ln  £i.  +  i  S  =  Ln  ^f  +  i  /FTh 
um     ^        Jh+1 


or 


Um  .  «  .       y/h+I 


-  T?  =       V"  '   —  (9) 


-v/f^i  +  Tf+h 


Before  proceeding  further,  It  should  be  pointed  out  that  the 
parameters  (f ,  h,  k)  could  have  been  selected  at  other  points  or  special 
use  could  have  been  made  of  the  symmetry  of  the  points  1  and  D.  A  care- 
ful analysis  of  the  various  possibilities  hat  shown  that  the  t-plane 
selected  herein  will  yield  the  least  amount  of  difficulty  In  the  evalua- 
tion of  the  integrals  to  be  performed  later. 

The  complex  potential  w  can  be  expressed  in  ten*  of  t  by  the  method 
of  sources  and  sinks  as 


U*c       U-a  U  •• 

W  ■  -£-  Lnt  -  -§-  Ln(t-k)  -  -«-—  Ln(t+f)         (10) 
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From  this  relationship  and  from  the  definition  of   C  ,  one  obtain; 


z  =  '^ti  JL^+1  +  v't_hJ  Lt  -  &  -  f?fJ dt         <1" 


The  integrals  appearing  in  the  above  equation  may  be  evaluated  for 
various  intervals  in  the  t-plane  or  for  various  distances  between  various 
points  in  the  physical  plane.   The  position  of  the  point  B  relative  to 
the  point  F  is  given  by: 

-t  +  ia  =  !  z  | 

t=h 
which,  upon  evaluation  in  a  straightforward  manner,  yields 

|  =    -  -4-  [-  |  Ln  4Tl    +   1    +    fiATl   Ln  ^1    +  fo1 
*vh+l  v^h+l    -    1  v'k+l    -  v'h-l 

(12) 

c    mKT       *-  e     r— ,     L_    -1     /fTl 


+  |  ^h"  -   2T7  J  /f-1    tan   x   ,^=i  -   irrj  §  /f+h 


It   should  be  noted   that  c/a,     6     ,    e/a,    and     T?  are   related   by   the   equa- 
tion of  continuity   as 

f  =    6   +  f   71  (13) 

The   evaluation   of   the    integrals    indicated  in   Eq.    (11)    from   t   ■   - i 

to   t   =   -  °°    and   the   separation   of   the   real   and  imaginary   parts   similarly 
yields 


f  =    -  -±—   I  §  n    -    * 6/k+l    +  £  V,T+h   Ln  ^  f+h    +   v 1+h 
3  <h+l    ua  a  v'f+h    -  ,'l+h 


c     c~  _      vl+h   +  vh         ~«     £ — r  -1        k-h 

-  7  vh  Ln —   -    26  vk-h    tan  t~4t 

a  'i  Tv,  V,  "    l+n 


(14) 
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It  is  thus  apparent  that  the  velocity  ratios  and  the  relative  setbacks 
t/a  and  s/a  and  the  angle   jS   may  be  evaluated  by  assigning  various 
values  to  the  parameters  f,  h,  and  k  in  the  t-plane.   It  is  further 
obvious  that  the  width  of  the  jet  at  C  may  be  evaluated  without  diffi- 
culty by  carrying  out  similar  integrations  from  D  to  C  and  from  B  to  C. 
Since  the  jet  width  or  the  contraction  coefficient  is  not  needed  in  the 
present  analysis,  the  results  concerning  the  contraction  coefficient  will 
not  be  presented  here. 

The  expressions  for  t/a  and  s/a  are  of  such  nature  that  one  cannot 
assign  special  values  to  them  and  expect  to  find  the  corresponding 
values  of  f,  g,  and  k  and  the  angle  3  .   The  nonlinearity  of  the  equa- 
tions makes  it  mandatory  that  one  assign  arbitrary  values  to  the  para- 
meters in  the  t-plane  in  a  fairly  systematic  manner  and  obtain  the 
corresponding  values  of  t/a,  s/a,  and  P   .    In  this  manner,  one  obtains 
an  infinite  variety  of  possible  solutions  and  then  selects  those  which 
conform  to  the  selected  geometry.   This  procedure  which  is  common  to  all 
free-streamline  problems,  has  been  applied  to  the  present  problem  through 
the  use  of  the  IBM  360  computer  and  the  results  obtained  for  both 
t/a=-l  ,  s/a  -  0,  e/a»  l  and  for  t/a  -  0.5  e/a  -  1,  and  s/a  ■  0  are 
plotted  in  terms  of  the  jet  deflection  angle  0!   versus  the  ratio  of 
the  velocities.   These  will  be  discussed  later  in  connection  with  the 
discussion  of  the  results.   The  information  for  other  combinations  of 
t/a,  s/a,  and  e/a  is  available  in  the  form  of  computer  print-out  and 
they  will  not  be  reproduced  here  for  it  is  practically  impossible  to 
cover  the  entire  range  in  the  form  of  a  simple  plot.   Thus  only  those 
cases  which  are  of  practical  significance  are  presented  herein  and 
compared  with  those  obtained  experimentally.   The  computer  program  which 
will  generate  solutions  for  all  ranges  of  the  parameters  encountered  in 
the  analysis  is  presented  in  the  Computer  Output  Section. 
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III.   EXPERIMENTAL  EQUIPMENT 

The  experimental  apparatus  consisted  of  the  }et  assemhlv,  supply 
system,  hot  wire  anemometer  system,  and  the  velocity  calibrator. 

Jet  Assembly 

The  jet  assembly  was  made  of  plexiglass  for  ease  of  machining  and 
flow  visualization.   As  shown  in  Fig.  2,  two  ten  inch  long,  one  inch 
deep  and  one-eighth  inch  wide  channels  were  cut  out  and  sandwiched  be- 
tween two  one  half  inch  plexiglass  plates.   The  axes  of  the  two  jets 
(referred  to  hereafter  as  the  main  jet  and  auxiliary  jet)  were  normal  to 
each  other  since  in  the  majority  of  the  momentum -exchange  amplifiers 
these  two  jets  are  made  to  intersect  each  other  at  an  angle  of  90 
degrees.   In  order  to  prevent  premature  wall-attachment  via  Coanda  ef- 
fect, the  outer  edges  of  the  two  jets  were  sharply  curved  back.   The 
three-eighths  inch  diameter  supply  lines  were  connected  to  the  inlets 
of  the  two  channels  and  a  smooth  transition  was  provided  between  the 
circular  and  the  rectangular  sections.   In  addition,  several  wire  screens 
were  placed  near  the  transition  region  for  the  purpose  of  eliminating 
secondary  swirling  flows  which  might  have  developed  in  the  long  circular 
supply  line.   The  outer  wall  of  the  auxiliary- jet  channel  was  made  mov- 
able for  the  purpose  of  adjusting  the  relative  set-back  (see  Fig.  2). 
All  channel  surfaces  were  polished  carefully  and  the  dimensions  were 
set  with  gage  blocks  to  an  accuracy  in  the  order  of  0.001". 

Supply  System 

The  working  fluid  (air)  was  supplied  from  a  50  Hp  compressor. 
Air  was  then  passed  through  a  microfilter  for  the  purpose  of  removing 
oil,  dust,  and  other  impurities  and  was  directed  to  pressure  regulators 
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and  Fischer-Porter  rotometers.   Both  the  main  jet  and  the  auxiliary  jet 
were  supplied  air  independently  from  the  main  flow  line  coming  out  of 
the  primary  pressure  regulator.   Each  jet  had,  in  addition,  a  series  of 
pressure  regulators  ranging  from  ten  to  130  psig  and  an  independent  flow- 
rator  whose  capacity  was  identical  with  that  of  the  other  flowrator.   The 
maximum  capacity  of  each  flowrator  was  14.6  SCFM  at  14.7  psia  and  70  de- 
grees with  a  3/4-27G-10/83  tube  and  3/4-GNSVGT-54  float.   Finally,  in 
order  to  determine  the  actual  flow  rate  passing  through  the  flowrators 
at  pressures  other  than  14.7  psia,  each  flowrator  was  equipped  with  a 
pressure  gage  and  the  flow  rates  were  obtained  by: 

Q  ■  flowrator  meter  reading  x  100%  full  flow  xCF(P) 
where 

CF(P)  is  the  pressure  correction  factor  and  is  provided  by  Fischer- 
Porter  Company  in  the  form  of  a  chart  for  each  flowrator. 
Hot-Wire  Anemometer  and  Calibrator 

The  Model  1050  constant  temperature  hot-wire  anemometer  of  the  Therwo- 
Systems,  Inc.,  was  used  in  the  measurement  of  velocities  end  turbulence 
intensities.   The  hot-wire  probe  holder  was  mounted  on  a  micrometer 
block  capable  of  providing  a  translation  (along  a  slot  along  the  edges 
of  the  top  and  bottom  cover  plates)  and  a  rotation  (in  the  mid-plane  of 
the  jets)  to  the  sensor  holder.   The  micrometer  enabled  the  measurement 
of  distances  along  any  section  in  the  mid-plane  between  the  top  and 
bottom  cover  plates  to  an  accuracy  of  1/160". 

The  probe  used  in  the  tests  had  a  cold  resistance  of  6.45  ohms  end 
an  operating  resistance  of  9.70  ohms.   The  probe  was  celibrated  with  the 
Model  1125  calibrator  of  Thermo-Systems,  Inc.   All  calibrations  were 
carried  out  with  the  high-range  calibration  orifice  of  the  calibrator 
which  had  a  0.150"  diameter. 
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In  order  to  calculate  the  velocity  of  the  flow  through  the  orifice, 
the  pressure  of  the  calibration  chamber  was  connected  to  water  mano- 
meters of  three  different  ranges  to  augment  the  accuracy  of  the  calibra- 
tion. 

A  schematic  of  the  experimental  layout  is  shown  in  Fig.  4.   Figures 
5  and  6  are  the  photographs  of  the  experimental  system,  the  calibrator, 
micrometer  assembly,  and  the  plexiglass  test  model. 
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Plexiglass  jet  assembly. 


Model  1125  calibrator. 


FIGURE  6.  JST  ASSEMBLY  AIJD  CALIBRATOR, 
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IV.   EXPERIMENTAL  PROCEDURE 

The  experiments  consisted  of  the  setting  of  proper  dimensions  to 
the  let  channels,  the  calibration  of  the  hot-wire  probe,  setting  of  the 
anemometer,  setting  of  the  selected  flow  rates  for  each  jet,  the  deter- 
mination of  the  let  deflection  angle  01      through  the  use  of  a  flow  visuali- 
zation method,  and  the  measurement  of  velocities  and  turbulence  intensities 
along  lines  normal  to  the  jet  axis  in  the  mid-plane  of  the  jet  at  various 
distances  from  the  nozzle  exit. 

Setting  of  the  Channel  Dimensions 

It  is  apparent  that  an  infinite  variety  of  geometrical  configura- 
tions may  be  tested  even  if  the  jet  impingement  angle  is  maintained  at 
90  degrees.   For  this  reason,  attention  was  given  only  to  those  geometries 
which  have  considerable  practical  significance  as  far  as  the  momentum- 
exchange  amplifiers  are  concerned.   In  the  majority  of  such  devices,  the 
channel  widths  are  kept  identical  and  the  side-wall  of  the  auxiliary  jet 
is  set  back  a  distance  equal  to  either  the  nozzle  width  or  one-half  of 
the  nozzle  width.   Thus,  in  the  present  study,  the  channel  widths  were 
kept  identical  (1/8")  and  two  different  set-backs  were  used  for  the  auxi- 
liary jet,  i.e.,  t/a  was  chosen  in  the  first  series  of  tests  equal  to 
0.5  and  in  the  second  series  of  tests  equal  to  unity.   All  dimensions 
were,  as  cited  before,  set  with  gage  blocks  to  an  accuracy  of  0.001". 
Calibration  of  the  Hot-Wire  Probe 

The  hot-wire  probe  was  placed  1/16"  above  the  axis  of  the  calibra- 
tion orifice  as  recommended  by  the  manufacturer  [5]  and  the  pressure 
tap  of  the  calibrator  was  connected  to  a  water  manometer.   In  order  to 
increase  the  accuracy  of  the  pressure  measurement,  three  different 
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manometers  were  used.   These  consisted  of  manometers  0  to  0.2"  range 
with  0.01"  increments,  0  to  2"  range  with  0.1"  increments,  and  a  mano- 
meter with  0-36"  range  with  0.1"  increments.   The  first  two  manometers 
were  of  inclined  type  for  further  increase  of  the  accuracy  in  the  reading 
of  the  water  elevations.   The  velocities  through  the  orifice  were  then 
calculated  through  the  use  of  the  following  equation: 


OR      V   p  v  ' 

where 

A?   is  the  manometer  readings  and   p   the  density  of  air.   The 
corresponding  voltage  output  on  the  hot-wire  anemometer  was  recorded 
for  each  manometer  setting  and  then  a  velocity  vs.  volts  "calibration 
curve"  was  plotted  as  shown  in  Fig.  7. 

The  calibration  curve  as  obtained  above  is  of  limited  use  because 
of  the  difficulty  it  introduces  in  the  evaluation  of  velocity  for  each 
anemometer  reading  corresponding  to  velocity  and  turbulence  measure- 
ments along  various  lines,  at  various  distances  from  the  nozzle,  and 
for  various  combinations  of  the  flow  rates.   It  is,  therefore,  obvious 
that  one  has  to  fit  one  or  more  polynomials  to  the  calibration  data 
and  calculate,  through  the  use  of  the  computer,  the  velocities  and 
turbulence  intensities  corresponding  to  each  anemometer  and  Root  Mean 
Square  meter  reading.   For  this  purpose,  various  degrees  of  polynomials 
were  tried  for  the  representation  of  the  entire  calibration  data  and  it 
was  found  that  two  4th  degree  polynomials  with  a  common  ordinate  and 
tangent  at  the  point  of  intersection  would  be  most  satisfactory.   The 
polynomials  chosen  were  of  the  following  type: 

_        _2       —  3       —  4 

U(E)  =  a1   +  a2AE  +  a3(AE)  +  a4(AE)  +  a5(AE)        (16) 
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where  AE  =  E  -  E   and  E   is  the  anemometer  output  at  zero  flow  velocity 
o      o 

and  U  is  the  temporal  mean  velocity  at  a  given  point.   The  coefficients 

a   of  the  above  polynomial  were  determined  through  the  use  of  an  IBM  360 

(see  Appendix  B) .   One  4th  degree  equation  with  coefficients  a  was  used 

for  differential  anemometer  output  equal  or  below  E  -  E  ■  1.460  and 

another  4th  degree  polynomial  with  coefficients  b  was  used  for  anemometer 

differential  outputs  greater  than  E  -  E  ■  1.460.   The  coefficients  a  and 

o  n 

b   are  tabulated  in  Appendix  B» 
n 

The  calibration  curve  described  above  was  checked  frequently  during 
experiments  (at  least  once  after  the  completion  of  the  velocity  and 
turbulence  measurements  for  a  given  flow  combination)  for  the  purpose  of 
ensuring  that  neither  the  probe  nor  the  anemometer  system  had  drifted 
from  their  original  conditions. 

Setting  of  the  Anemometer 

Prior  to  each  experiment,  the  anemometer  was  turned  on  and  set  to 
warm  up  for  about  one-half  hour.   Then  the  cold  resistance  of  the  hot- 
wire probe  was  checked  to  make  sure  that  it  had  not  deviated  from  its 
previously  measured  value.   Then  the  operating  resistance  was  set  in 
accordance  with  the  regular  operational  instructions  of  the  anemometer 
[5]  and  the  system  was  set  to  'RUN'  position.   One  of  the  output  channels 
of  the  anemometer  was  connected  to  an  oscilloscope  and  the  flow  was  turned 
on  and  increased  gradually  until  an  oscillation,  due  to  feedback  in- 
stability in  the  amplifier  circuitries,  was  encountered.   Then  the  'TRIM' 
on  the  anemometer  was  adjusted  and  the  oscillation  wee  eliminated.   This 
procedure  was,  of  course,  carried  out  for  all  velocities  encountered  in 
the  experiments  by  increasing  and  decreasing  the  flow  rate. 
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The  oscilloscope  was,  of  course,  observed  at  all  times  during  the 
actual  measurement  to  make  sure  that  no  system-generated  noise  inter- 
fered with  the  actual  turbulence  of  the  flow. 

Setting  of  the  Selected  Flow  Rates 

The  flow  rates  through  each  channel  were  selected  on  the  basis  of 
the  Reynolds  numbers  encountered  in  the  application  of  large  momentum- 
exchange  devices.   The  velocity  of  the  main  jet  was  maintained  at  173.200 
ft/sec  which  corresponded  to  a  flowrator  reading  of  607,  of  the  full  capa- 
city and  to  a  Reynolds  number  of  11300.   The  flow  rates  of  the  auxiliary 
iet  ranged  from  60%  to  107,  at  increments  of  10%.   These  flow  rates  and 
the  corresponding  velocities  and  Reynolds  numbers  are  tabulated  below. 

Flow  Rate  Velocity  ft/sec         Reynold  Number 

60%  173.200  11300 

507o  146.500  9550 

40%  113.200  7  380 

30%  84.096  5485 

20%  56.064  3660 

10%  28.032  1832 

During  the  experiments,  a  flow  rate  was  selected  for  the  auxiliary 
let  from  the  above  table  and  the  flow  through  each  iet  was  maintained 
constant  (within  the  limits  of  the  experimental  uncertainty)  with  the 
use  of  separate  pressure  regulators  described  previously. 
The  Determination  of  the  Jet  Deflection  Angle 

The  iet  deflection  angle  Of     was  determined  through  the  use  of 
several  different  procedures.   The  first  two  of  these  procedures  con- 
sisted of  the  use  of  lamp-black  and  oil  mixture  and  the  use  of  powdered 
color  paint. 
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These  gave  the  jet  deflection  angle  to  an  accuracy  within  +2°. 
The  third  method  which  yielded  the  jet  deflection  angle  to  a  finer 
degree  of  accuracy  was  the  use  of  the  measured  velocity  profiles. 

The  lamp-black  and  oil  were  mixed  at  more  or  less  arbitrary  quanti- 
ties and  the  inner  face  of  the  base  plate  was  carefully  painted  with  a 
brush.   Prior  to  painting,  the  valves  and  the  pressure  regulators  were 
adjusted  in  such  a  manner  that  when  the  main  flow  was  turned  on,  the 
flow  combination  would  correspond  to  the  desired  flow  rates  in  each 
channel.   Obviously  this  latter  procedure  was  easily  carried  out  by 
adjusting  the  flows  prior  to  the  painting  and  then  shutting  off  the 
main  valve  rather  than  the  secondary  valves.   Then  when  the  main  valve 
was  turned  on,  the  resultant  jet  partly  washed  the  surface  of  the  paint 
and  left  a  clear  print  of  the  jet  configuration.   It  is,  of  course, 
fully  realized  that  the  boundary  layer  may  yield  a  flow  configuration 
on  the  paint  somewhat  different  from  that  which  would  be  encountered  in 
the  mid-plane.   Nevertheless,  the  method  yielded  the  jet  deflection 
angles  which  were  fairly  satisfactory  as  far  as  the  preliminary  deter- 
mination of  the  velocity  and  turbulence  traversing  planes  was  concerned. 
The  jet  deflection  angle  was  measured  relative  to  the  main  jet  by  plac- 
ing two  graph  papers  above  and  below  the  test  section  and  by  lining  up 
the  paint  trace  with  the  lines  on  the  graph  papers. 

The  second  method  used  in  the  determination  of  the  jet  deflection 
angle  consisted  of  the  preparation  of  a  white  oil-painted  paper  and 
the  introduction,  into  the  supply  lines,  of  small  amounts  of  powdered, 
colored,  dry  paint  (different  color  for  each  channel).   Tha  painted  paper 
was  glued  to  the  inner  face  of  the  bottom  plate  with  rubber  cement.   The 
powdered  paint  released  from  the  supply  line  followed  the  Jet  configura- 
tion and  stuck  on  the  freshly  painted  paper.   Thus,  this  method  quickly 
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yielded  the  jet  configuration  and  the  deflection  angle  and  also  gave  a 
pretty  good  idea  about  the  mix-up  of  the  jets  in  the  interaction  region. 
In  all  the  experiments  cited  herein,  the  mixing  region  and  of  course  the 
remainder  of  the  jet  came  out  at  a  color  which  one  would  have  obtained 
by  mixing  them  thoroughly  in  a  separate  container.   A  sample  jet  con- 
figuration as  obtained  in  this  manner  is  shown  in  Fig.  8. 

The  Measurement  of  Velocity  and  Turbulence  Intensities 
The  velocities  were  measured  along  the  lines  normal  to  the  jet 
centerline  determined  through  one  of  the  flow  visualization  techniques 
cited  above.   The  hot-wire  probe  was  moved  at  suitable  increments  along 
lines  normal  to  the  jet  centerline  and  both  the  velocity  and  Root  Mean 
Square  (RMS)  meters  were  read  and  recorded.   The  velocities  and,  of 
course,  turbulence  intensities  were  measured  at  six  different  stations 
for  six  different  flow  combinations.   These  stations  were  located  at  0.15", 
0.3",  0.8",  1.5",  2.5",  and  4"  from  the  point  where  the  inner  walls  of 
the  two  jets  intersected.   As  will  be  mentioned  later,  only  34  velocity 
and  34  turbulence  intensity  traces  are  presented  instead  of  36  of  each. 
(Because  the  flow  rates  of  1/6  and  2/6  made  it  impossible  to  maneuver  the 
hot-wire  probe  at  the  first  station  (0.15"  from  the  nozzle)  due  to  the 
fact  that  the  traversing  line  fell  within  the  auxiliary  jet  channel.) 
Figure  9  is  a  schematic  drawing  of  the  velocity  traversing  line  rela- 
tive to  the  jet  channels. 


The  turbulence  intensity     V  u  /fj    was  calculated  by 

,_2     _ 

,/-  2      Me     dU 

;_2  - 

where  Ve    is  the  output  of  the  RMS  meter  and  dU/dE     is  the  deriva- 
tive of  Eq.  16  or,  in  other  words,  the  complement  of  the  slope  of  the 
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Ua/Um    (173.2/173.2)    =  1.000 
Re  (auxiojet        )    s=   1 1 300 

Ho  (main  jet        )   ==  11300 


Main  jet 


FIGURE  6.  VISUAL  JET  CONFIGURATION. 
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1 •Main  jot  channel 
2.Auxi.jet    channel 
3. Deflection  angle 
4-5  Sweeping  path 
6, Jet  centerline 
7. Station  distance (d) 


FIGURE  9.      SCHEMATIC  OF  VELOCITY  TRAVERSING  LINE 
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calibration  curve.   The  derivation  of  the  above  equation  is  shown  in 
Appendix  C.   The  turbulence  intensity  profiles  were  plotted  like  the 
velocity  profiles  by  means  of  an  IBM  360  computer  and  will  be  discussed 
subsequently. 
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V.   DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

The  purpose  of  the  present  investigation  was  the  understanding  of 
the  characteristics  of  the  impingement  of  two  plane  jets.   These  character 
istics  essentially  were  the  deflection  angle  of  the  resulting  iet,  and  the 
velocity  and  turbulence  distributions. 

The  deflection  angles  obtained  experimentally  are  shown  in  Fig.  10 
for  a  1/16"  set-back  (t/a  =  -0.5).   The  similar  data  for  a  1/8"  set-back 
(t/a  =  -1.0)  are  shown  in  Fig.  11.   Also  shown  on  both  figures  are  the 
results  obtained  theoretically  through  the  use  of  Eqs .  (12)  and  (14). 

It  is  apparent  from  the  examination  of  Figs.  10  and  11  that  the 
agreement  between  the  theoretical  and  experimental  results  is  quite 
satisfactory,  at  least  for  the  flow  rates  and  set-backs  used. 

The  velocity  distributions  for  various  angles  of  deflection  at 

various  distances  along  the  nozzle  exit  are  shown  in  Figs.  12-29  in  a 

normalized  form,  i.e.,   U/U  versus  x/a.   The  following  facts  become  ap- 

m 

parent  upon  a  brief  examination  of  these  figures: 

(a)  The  velocity  profiles  become  flatter  and  perfectly  sym- 
metrical as  the  distance  increases  from  the  nozzle  exit  to  the  section 
where  the  velocity  is  measured.   This,  of  course,  is  an  observation 
which  is  true  for  all  iet  flows. 

(b)  The  mixing  of  the  two  jets  takes  place  more  quickly  for 

the  smallest  iet  velocity  ratios  such  as  U  /U   =  1/6.   As  the  velocity 

a  m 

ratio  increases,  the  velocity  profile  exhibits  two  maximums  with  a 
trough  in  the  middle.   This  points  out  that  the  mixing  of  high  velocitv 
iets  requires  larger  distances  along  the  iet  but  certainly  a  distance 
less  that  approximately  d/a  =  6.   The  significance  of  this  observation 
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is  in  relation  with  the  theoretical  analysis.   As  has  been  cited  in 
the  development  of  the  theoretical  analysis,  the  free-streamline  velo- 
cities on  the  two  sides  of  the  inviscid  jet  had  to  be  assumed  to  be 
identical.   (The  velocity  profiles  for  all  velocity  ratios,  near  the 
nozzle  where  the  free-streamline  analysis  is  supposed  to  hold  true,  show 
an  almost  uniform  distribution  with  the  exception  of  two  peaks  and  a 
trough  whose  amplitude  is  about  1/8  of  the  maximum  of  the  normalized 
velocity.   The  velocity  profiles  for  each  angle  are  combined  in  Figs. 
30-35  for  the  purpose  of  giving  a  quick  idea  about  the  evolution  of  the 
let  from  the  nozzle  downstream. 

Figures  30-41  are  plotted  using  the  parameters  U/U  vs.  x/d.   The 
purpose  of  these  plots  was  to  explore  the  similarity  in  the  velocity 
profiles  in  the  initial  region  of  jet  development.   Previous  studies  [6] 
on  the  evaluation  of  single  two-dimensional  and  axi-symmetric  jets  have 
shown  that  plots  prepared  in  the  manner  cited  above  resulted  in  a  single 
normalized  curve.   This  was,  to  be  sure,  no  mixing  occurs  with  the  ex- 
ception of  that  with  the  surroundings.   The  plots  presented  herein  show 
that  such  a  similarity  exists  only  for  distance  of  d/a  >   2.4.   This  is 
evidenced  by  the  fact  that  the  two  velocity  profiles  at  distances  of 
d/a  ■  1.2  and  d/a  ■  2.4  do  not  fall  on  the  rest  of  the  velocity  profiles 
exhibiting  similarity.   This  may  be  interpreted  as  a  consequence  of  the 
mixing  of  two  jets  of  unequal  velocity. 

The  turbulence  intensities  were  plotted  as  a  function  of  x/a  for 
various  angles  of  delfection  and  for  various  stations  along  the  jet  and 
are  shown  in  Figs.  42-59.   The  main  characteristic  of  all  the  turbulence 
distributions  is  that  there  are  two  high  regions  of  turbulence  and  a 
region  of  relatively  low  turbulence  in  the  central  core  of  the  jet.   The 
result  is  not  unexpected  in  view  of  the  fact  that  •«ch  Jet  has  the 
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highest  shear  and  mixing  with  the  surrounding  ambient  flow  at  its  two 
sides.   It  is  also  apparent  from  these  figures  that  the  intensity  of 
turbulence  increases  with  increasing  distances  along  the  jet.   This  is 
not  so  much  because  of  the  absolute  increase  of  the  RMS  value  of  the 
local  values  of  velocity  fluctuations,  but  rather  because  of  the  rela- 
tive decrease  of  the  magnitude  of  velocity  along  the  jet.   The  turbulance 
intensity  profiles  for  each  angle  are  combined  in  Figs.  60-65. 
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VI.   RECOMMENDATIONS  FOR  FURTHER  WORK 

The  following  are  recommendations  for  further  studies: 

a.  Jet  deflection,  velocity  and  turbulance  distributions  be 
studied  for  other  set-back  and  channel  widths. 

b.  The  study  of  the  above  parameters  be  extended  to  lower 
Reynolds  numbers  where  one  or  both  of  the  jets  are  laminar. 

c.  A   theoretical  analysis  be  carried  out  for  unidentical 
free-stream  velocities,  and  finally, 

d.  The  aspect  ratio  (the  ratio  of  the  width  of  the  channel 
to  the  depth  of  the  channel)  be  changed  from  2  to  8. 
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APPENDIX  A 

Evaluation  of  the  Integrals 
The  integrals  indicated  by  Eq.  11  may  be  evaluated  as  follows 


=    '  c^+l  dt  =  c   {  2ytTi  .  Ln  ^M.  +71  } 


I2    =    -    a6 


/y/t+1 

t-k 


dt   = 


-   a6    {   2/bTl    -  Vk+I   Ln  ^   +  ^ 

Jt+i  -  <s/k~+i 


I3   -  -  el) 


t+1 
t+f 


dt 


=   -   eT)  {   2/t+l   +   2-v/f^l    tan' 


-1    /f-1 
t+1 


J4    "   C 


^^  dt   =  c    {   2^/t^h   +   2^   tan"1   J-£g  J 


T-h 


I5   =    -   a6        ^^  dt   =   -aft    (2/t^h   -  Jk-h 


Ln 


i6  =  -  erj 


*g=k  dt   =   -el?  {2/t^h  +   2^/f+h    tan"1   /"f^-  } 
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The  above  results  are  evaluated  between  proper  limits  to  yield  the 
normalized  distances  s/a  and  t/a  given  by  Eqs.  (12)  and  (14).   Note 
that  regardless  of  the  integration  limits  the  sum  of  the  first  terms 
in  the  integrated  expressions  is  zero  because  of  continuity  (see  Eq. 
(13)). 
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APPENDIX  B 

Determination  of  the  Coefficients  of  the  Polynomial 
One  of  the  computer  programs  for  the  determination  of  the  coeffici- 
ents of  the  polynomials  describing  the  two  parts  of  the  calibration  curve 
is  shown  on  page 

The  polynomials  for  either  part  of  the  calibration  curve  were  as- 
sumed to  be  of  the  following  type. 

U(E)  =  a1   +  a2£E  +  a3(AE)2  +  a4  (Ae)  3  +  a5(AE)4 

The  coefficients  a  and  b   for  the  two  polynomials  are  tabulated  below: 

n      n  J 

b 

n 

bx  -  3728.036 

b2  —9160.219 

b3  -  8384.484 

b,  —3353.178 
4 

b5  -  512.192 

The  second  program  used  in  the  determination  of  the  coefficients  is 
not  presented  herein  for  it  was  identical  in  every  aspect  with  the  pro- 
gram used  for  the  first  part  of  the  calibration  curve  with  the  excep- 
tion of  course,  of  the  test  data  and  the  range  of  iteration.   The 
resulting  calibration  curve  and  the  test  data  are  shown  in  Fig.  7. 


a 

n 

al 

=  0.764 

a2 

=  17.491 

a3 

-  20.558 

a4 

=-20.184 

a5 

-  24.969 
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APPENDIX  C 

Analysis  of  the  Turbulence  Measurement 

Evaluation  of  the  intensity  of  the  turbulence  is  based  on  the  well- 

2  2       ^M 

known  King  s  Law  given  as  E  -E   ■  BU   where 

E   is  the  anemometer  output  at  zero  velocity,  E  is  the  output 
o 

for  the  velocity  U,  and  B  is  a  constant.   The  exponent  n  is  ordinarily 
assumed  to  be  2.   In  most  of  the  hot-wire  work  covering  a  large  range 
of  velocity,  the  exponent  may  deviate  considerably  from  2  and  it  is 
for  this  reason  that  in  precision  work  one  must  use  a  calibration  curve. 
The  evaluation  of  such  a  calibration  curve  was  described  in  the  pre- 
vious appendix.   In  the  following  the  evaluation  of  the  corresponding 
turbulence  intensities  will  be  described. 

For  a  turbulent  flow  the  instantaneous  value  of  the  anemometer  out* 

-    f^ 
put  is  given  by  E  =  E  +  Ve 


Where   E   is  the  voltage  output  corresponding  to  the  mean  velocity  at 

/T2 
a  given  point  and  M   e      is  the  Root  Mean  Square  value  of  the 

IZ1 

voltage  output  corresponding  to  V  u      (RMS  value  of  the  fluctuating 
component  of  velocity). 
Writing 

E   =   f  (U)  ,         dE   =    f '  (U)    dU 

and 

U   =  U   +Vu 


or 


Z2_    ._.      in  ,  _. 


/    2  _  '     2  _ 

f  (U)     =    f  (U    +  Ju     )  =    f  (U)     +  7u     f    (U) 
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and  noting  that 

-2        J? 


u 


f  (U) 
Thus  the  intensity  of  turbulence  reduces  to 


_2        -2       /-2   - 
u       Ve       Ve   dU 


U      Uf  (U)      U    dE 

where 

f'(U)  =   € 
dU 

CI 
The  turbulence  intensity   ^u  y-r=    was  evaluated  once  again  through 

Hi 
the  use  of  computer  by  feeding  in  the  measured  values  of   ^e     ,  U, 

and  the  values  of    ^VdE      as  ODtaine<^  from  the  two  polynomials  for 
respective  parts  of  the  calibration  curve. 
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